Improvement of Distance Measurement Based on Dispersive Interferometry Using Femtosecond Optical Frequency Comb

Since the dispersive interferometry (DPI) based on optical frequency combs (OFCs) was proposed, it has been widely used in absolute distance measurements with long-distance and high precision. However, it has a serious problem for the traditional DPI based on the mode-locked OFC. The error of measurements caused by using the fast Fourier transform (FFT) algorithm to process signals cannot be overcome, which is due to the non-uniform sampling intervals in the frequency domain of spectrometers. Therefore, in this paper, we propose a new mathematical model with a simple form of OFC to simulate and analyze various properties of the OFC and the principle of DPI. Moreover, we carry out an experimental verification, in which we adopt the Lomb–Scargle algorithm to improve the accuracy of measurements of DPI. The results show that the Lomb–Scargle algorithm can effectively reduce the error caused by the resolution, and the error of absolute distance measurement is less than 12 μm in the distance of 70 m based on the mode-locked OFC.


Introduction
With the development of science and technology, distance measurement with higher precision is widely used in modern industry and scientific research, such as industrial production, measurement of gravitational waves, acquisition of black-hole images, and observation of satellite formations. However, due to the current technical level, it is very challenging to realize distance measurement with a large length and high accuracy at the same time, so it is urgent to find a method with high-precision and long-distance range.
The emergence of OFCs has opened up a new way to solve this problem. The related concepts were first proposed by Hänsch in 1978 [1]. OFCs are a kind of pulse-light that have a time width in femtoseconds [2]. They provide a new research direction for so many fields of distance measurement, precision spectroscopy frontier physics, and so on [3,4]. They have been closely related to the development of the country and people's lives.
In 2000, Minoshima and Matsumoto of the National Institute of Metrology in Japan applied an OFC for absolute distance measurement for the first time [5]. This kind of distance measurement method based on OFCs includes synthetic wavelength interferometry [6][7][8], multi-wavelength interferometry [9], time-of-flight [10], dispersive interferometry [11], and optical sampling method [12].
In 2006, Joo and Kim et al. proposed a scheme for absolute distance measurement by DPI, and achieved a 1.46 mm non-ambiguity range and a resolution of 7 nm in 0.89 m [13,14], which also ignored the error brought by the non-uniform frequency sampling of the spectrometer. In 2013, Shujian Xing used a combination of multi-wavelength interferometry and time-of-flight based on a Michelson interferometer to conduct distance measurement , sin(π f r t) .
By using the above two equations, we can obtain E(t) from Equation (8) as follows: where p = π f r t. The modulated signal can be simplified as follows: E m (t) = sin[(M + 1)p] sin(p) . (12) By using Equation (12), the expression of the OFC E(t) in Equation (8) can be rewritten as Equation (7), which indicates that the value of M cannot influence the form of E(t), and Equation (7) can be used as the mathematical model of the OFC.

Characteristics of the Frequency Comb
By analyzing the expression (11) and (12), we can know that E(t) is composed of a modulated signal and a carrier signal. Moreover, the modulated signal E m (t) determines the characteristics of E(t). The intensity of E 1 (t) is equal to the intensity of E m (t), and we can obtain its peak intensity by calculating the limit of its modulated signal when p = πk (k is an integer number).
Ideally, the peak intensity of the OFC depends on the number of its frequency components only. The modulated signal is periodic, and we can get its period by calculating the distance of two neighbouring crests. From Equation (13), the period of the modulated signal is ∆p = πf r ·∆t = π, and the period T can be derived as T = 1/f r . The width of the pulse in the time domain can be derived from Equation (6). When sin[(M + 1)πf r t] = 0 and sin(πf r t) = 0, there are two closest zero-crossing points z 1 and z 2 on both sides of the crest. The distance of these two points can be calculated as follows: where W can be used to estimate the width of the OFC in the time domain. As shown in Figure 1, we can obtain the simulated signal based on Equation (7). Here, the value of M is 381, f r = 10 GHz, and f c = 194.2 THz. The amplitude of the pulse peak is 382 and the period of the pulse train is 0.1 × 10 −9 s, which meets the conclusions of Equation (13). Moreover, the single pulse circled by the red box is also shown in Figure 1b and, according to Equation (14), the full width at half maximum of the pulse in the time domain is about 5.24 × 10 −13 s.
Ideally, the peak intensity of the OFC depends on the number of its frequency components only. The modulated signal is periodic, and we can get its period by calculating the distance of two neighbouring crests. From Equation (13), the period of the modulated signal is ∆p = πfr•∆t = π, and the period T can be derived as T = 1/fr. The width of the pulse in the time domain can be derived from Equation (6). When sin[(M + 1)πfrt] = 0 and sin(πfrt) ≠ 0, there are two closest zero-crossing points z1 and z2 on both sides of the crest. The distance of these two points can be calculated as follows: where W can be used to estimate the width of the OFC in the time domain. As shown in Figure 1, we can obtain the simulated signal based on Equation (7). Here, the value of M is 381, fr = 10 GHz, and fc = 194.2 THz. The amplitude of the pulse peak is 382 and the period of the pulse train is 0.1 × 10 −9 s, which meets the conclusions of Equation (13). Moreover, the single pulse circled by the red box is also shown in Figure 1b and, according to Equation (14), the full width at half maximum of the pulse in the time domain is about 5.24 × 10 −13 s.

The Principle of Dispersive Interferometry
The intensity of the interference signal can be written as where Erp(t) and Emp(t) denote the Fourier transform of the reference pulse train (Erp(t)) and the measurement pulse train (Emp(t)), respectively. Here, we assume that the measurement signal and reference signal have the same amplitude, the measurement signal can be considered as the reference signal with a time delay τ, and it can be written as Emp(t) = Erp(t−τ). According to the time shift property of Fourier transform, the relation of these two signals can be shown as follows: Furthermore, the interference signal I(f) can be simplified as 2 2 ( ) 2 ( ) 2 ( ) cos (2 ) rp rp In the frequency domain, I(f) is modulated by a signal whose frequency is τ. The distance difference L between the reference signal and the measurement signal can be calculated as L = τc/2.  The intensity of the interference signal can be written as where E rp (t) and E mp (t) denote the Fourier transform of the reference pulse train (E rp (t)) and the measurement pulse train (E mp (t)), respectively. Here, we assume that the measurement signal and reference signal have the same amplitude, the measurement signal can be considered as the reference signal with a time delay τ, and it can be written as E mp (t) = E rp (t − τ). According to the time shift property of Fourier transform, the relation of these two signals can be shown as follows: Furthermore, the interference signal I(f ) can be simplified as In the frequency domain, I(f ) is modulated by a signal whose frequency is τ. The distance difference L between the reference signal and the measurement signal can be calculated as L = τc/2.
Based on the conclusion of Equation (7), we obtain two simulated signals in the time domain, E rp (t) and E mp (t). To avoid the distortion of the modulated signal which is extracted by selecting the peak value of each frequency component of the interference signal, we set E mp (t) = 2E rp (t − τ), so that the modulation depth of the interference signal can be shallower and the peak values are above zero. The interference signal I(f ) is shown in Figure 2. The values of the center frequency f c , repetition frequency f r , and M used in our program are 194.2 THz, 10 GHz, and 381, respectively. The time difference between the reference comb and the measurement comb is 4 × 10 −12 s, as shown in Figure 3b, and the frequency of the modulated signal is 4 × 10 −12 Hz, which is equal to the value of the time delay.
Based on the conclusion of Equation (7), we obtain two simulated signals in the time domain, Erp(t) and Emp(t). To avoid the distortion of the modulated signal which is extracted by selecting the peak value of each frequency component of the interference signal, we set Emp(t) = 2Erp(t − τ), so that the modulation depth of the interference signal can be shallower and the peak values are above zero. The interference signal I(f) is shown in Figure 2. The values of the center frequency fc, repetition frequency fr, and M used in our program are 194.2 THz, 10 GHz, and 381, respectively. The time difference between the reference comb and the measurement comb is 4 × 10 −12 s, as shown in Figure 3b, and the frequency of the modulated signal is 4 × 10 −12 Hz, which is equal to the value of the time delay.

The Resolution of Dispersive Interferometry
The absolute distance measurement with dispersive interferometry depends on the frequency resolution of the interference signal. According to the Nyquist condition for sampling, the resolution of the time delay can be shown as follows: where B is the spectral width of the frequency comb, Fs is the sampling frequency, and ∆t is the sampling interval. Based on equality L = τc/2, the resolution of distance can be shown as follows: Based on the conclusion of Equation (7), we obtain two simulated signals in the time domain, Erp(t) and Emp(t). To avoid the distortion of the modulated signal which is extracted by selecting the peak value of each frequency component of the interference signal, we set Emp(t) = 2Erp(t − τ), so that the modulation depth of the interference signal can be shallower and the peak values are above zero. The interference signal I(f) is shown in Figure 2. The values of the center frequency fc, repetition frequency fr, and M used in our program are 194.2 THz, 10 GHz, and 381, respectively. The time difference between the reference comb and the measurement comb is 4 × 10 −12 s, as shown in Figure 3b, and the frequency of the modulated signal is 4 × 10 −12 Hz, which is equal to the value of the time delay.

The Resolution of Dispersive Interferometry
The absolute distance measurement with dispersive interferometry depends on the frequency resolution of the interference signal. According to the Nyquist condition for sampling, the resolution of the time delay can be shown as follows: where B is the spectral width of the frequency comb, Fs is the sampling frequency, and ∆t is the sampling interval. Based on equality L = τc/2, the resolution of distance can be shown as follows:

The Resolution of Dispersive Interferometry
The absolute distance measurement with dispersive interferometry depends on the frequency resolution of the interference signal. According to the Nyquist condition for sampling, the resolution of the time delay can be shown as follows: where B is the spectral width of the frequency comb, F s is the sampling frequency, and ∆t is the sampling interval. Based on equality L = τc/2, the resolution of distance can be shown as follows: Thus, we can conclude that the spectral width of the frequency comb decides the resolution of the dispersive interferometry: the wider spectrum of the frequency comb, the higher the resolution. The simulation of measurement resolution is shown in Figure 4. The values of the center frequency f c , repetition frequency f r , and M used in our program are 194.2 THz, 10 GHz, and 381, respectively, and thus the spectral width of the OFC is 3.8 THz. According to Equation (19), we can obtain the resolution of distance as 0.078 mm, which conforms to the simulation results.
Thus, we can conclude that the spectral width of the frequency comb decides the resolution of the dispersive interferometry: the wider spectrum of the frequency comb, the higher the resolution.
The simulation of measurement resolution is shown in Figure 4. The values of the center frequency fc, repetition frequency fr, and M used in our program are 194.2 THz, 10 GHz, and 381, respectively, and thus the spectral width of the OFC is 3.8 THz. According to Equation (19), we can obtain the resolution of distance as 0.078 mm, which conforms to the simulation results.

Non-Ambiguity Range of Dispersive Interferometry
The sampling interval of the interference signal and each frequency width of the optical frequency comb decides the measurement range of distance. The spectral width of the frequency comb and the resolution of the spectrometer also decide the sampling interval. When the spectrometer can distinguish each frequency of the frequency comb, the frequency comb can be seen as a sampled signal with an equal sampling interval in the frequency domain. For restoring the sampled signal correctly, the sampling frequency must be two times higher than the highest frequency in the sampled signal. The frequency comb can be seen as a signal with the sampling frequency 1/fr, so the highest frequency of the restored signal is 1/2fr. The non-ambiguity range is shown as follows: where LNAR is the non-ambiguity range and Lpp is the optical path difference between two adjacent pulses of the OFC. Because the width of each frequency mode of OFC affects the coherence of the measurement light and the reference light, the coherence length Lc determined by the linewidth of each frequency can be expressed by the following equation: The ∆f is the linewidth of each frequency. To prevent a measurement dead-zone we must meet Lc > LNAR. Thus, we can adjust the size of fr to realize the absolute distance measurement without a dead zone. The simulation of the non-ambiguity range is shown in Figure 5:

Non-Ambiguity Range of Dispersive Interferometry
The sampling interval of the interference signal and each frequency width of the optical frequency comb decides the measurement range of distance. The spectral width of the frequency comb and the resolution of the spectrometer also decide the sampling interval. When the spectrometer can distinguish each frequency of the frequency comb, the frequency comb can be seen as a sampled signal with an equal sampling interval in the frequency domain. For restoring the sampled signal correctly, the sampling frequency must be two times higher than the highest frequency in the sampled signal. The frequency comb can be seen as a signal with the sampling frequency 1/f r , so the highest frequency of the restored signal is 1/2f r . The non-ambiguity range is shown as follows: where L NAR is the non-ambiguity range and L pp is the optical path difference between two adjacent pulses of the OFC. Because the width of each frequency mode of OFC affects the coherence of the measurement light and the reference light, the coherence length L c determined by the linewidth of each frequency can be expressed by the following equation: The ∆f is the linewidth of each frequency. To prevent a measurement dead-zone we must meet L c > L NAR . Thus, we can adjust the size of f r to realize the absolute distance measurement without a dead zone. The simulation of the non-ambiguity range is shown in Figure 5: According to Equation (20), the time delay is 15 mm, which is equal to the resu the simulation. We can assume the measurement result is L. When no frequency mi occurs, such as in point A, the real value L1 is equal to L. However, in point B the m urement result is L2 and frequency mixing occurs, so the real value is in point D. In p C the measurement result is L3 and frequency mixing also occurs, so the real value point E. We can use the following equations to obtain the real value: where l is the optical path difference between the measuring signal and the reference nal, N is an integer which is obtained by rounding down the quotient of the distance t measured and Lpp/2, and Lx is the calculation results.
In actual calculation, the derivative of the curve near point B is less than 0. When actual distance increases, the measurement results decrease. The derivative of the c at point C is greater than 0. When the actual distance increases, the measurement re increase. Therefore, Equation (22) or Equation (23) can be selected according to the m urement results at a certain point and the measurement results nearby.
Since the light travels back and forth in the optical path, the actual measurem length is magnified to 2 times. For the convenience of expression, the resolution and ambiguity range described below refer to the measurement length of 1 time. Figure 6 shows the layout of the dispersive interferometry experiment. The source was a One-Five Origami-15 and the values of the center frequency fc and repet frequency fr used in our program were 192.3 THz and 250 MHz, respectively, while spectral width of the frequency comb was 15 THz. The spectrometer was a YOKOGA AQ6370D. The distance resolution of the system was about 10 μm, and the non-ambig range was about 0.6 m. According to Equation (20), the time delay is 15 mm, which is equal to the result of the simulation. We can assume the measurement result is L. When no frequency mixing occurs, such as in point A, the real value L 1 is equal to L. However, in point B the measurement result is L 2 and frequency mixing occurs, so the real value is in point D. In point C the measurement result is L 3 and frequency mixing also occurs, so the real value is in point E. We can use the following equations to obtain the real value:

Experiment with the Mode-Locked Frequency Comb
where l is the optical path difference between the measuring signal and the reference signal, N is an integer which is obtained by rounding down the quotient of the distance to be measured and L pp /2, and L x is the calculation results. In actual calculation, the derivative of the curve near point B is less than 0. When the actual distance increases, the measurement results decrease. The derivative of the curve at point C is greater than 0. When the actual distance increases, the measurement results increase. Therefore, Equation (22) or Equation (23) can be selected according to the measurement results at a certain point and the measurement results nearby.
Since the light travels back and forth in the optical path, the actual measurement length is magnified to 2 times. For the convenience of expression, the resolution and nonambiguity range described below refer to the measurement length of 1 time. The long distance absolute measurement experiment was carried out on an 8 derground optical guide rail of the National Institute of Metrology, China. The tes ror was placed on the air floatation vehicle and moved steadily on the guide rai experiment, a number of different measuring points were selected for measuremen a long distance, and the measurement result of the interferometer was used as th ard value. At the same time, the field environmental parameters were measured measurement process, the ambient temperature was 21.20 °C, the pressure was hPa, the air humidity was 25.17%, and the air refractive index ng obtained by th equation was 1.00026799.

Experiment with the Mode-Locked Frequency Comb
In the experiment, the laser interferometer (Agilent 5519b) was set to zero at tial position, and the frequency-domain interference signal of the frequency co recorded. The corresponding distance was recorded as L0. During the experim measured mirror was moved continuously, several measuring points were sele measurement, and the interferometer's signal and the interference signal of the fr comb were recorded at the same time. If the result of the interference signal of quency comb is LD, the optical path difference between the measured position starting position can be expressed as

Lomb-Scargle Algorithm and Experimental Results
The interference spectra at different positions were obtained in the measurem cess as shown in Figure 7. The data in the upper left corner are the distance betw measurement signal and the reference signal, and the data in the upper right co the measurement result obtained by the interferometer. The larger the interval b the two signals, the higher the carrier frequency of the interference signal and the the interference fringes. The long distance absolute measurement experiment was carried out on an 80 m underground optical guide rail of the National Institute of Metrology, China. The tested mirror was placed on the air floatation vehicle and moved steadily on the guide rail. In this experiment, a number of different measuring points were selected for measurement across a long distance, and the measurement result of the interferometer was used as the standard value. At the same time, the field environmental parameters were measured. In the measurement process, the ambient temperature was 21.20 • C, the pressure was 1016.10 hPa, the air humidity was 25.17%, and the air refractive index n g obtained by the Edlen equation was 1.00026799.
In the experiment, the laser interferometer (Agilent 5519b) was set to zero at the initial position, and the frequency-domain interference signal of the frequency comb was recorded. The corresponding distance was recorded as L 0 . During the experiment, the measured mirror was moved continuously, several measuring points were selected for measurement, and the interferometer's signal and the interference signal of the frequency comb were recorded at the same time. If the result of the interference signal of the frequency comb is L D , the optical path difference between the measured position and the starting position can be expressed as

Lomb-Scargle Algorithm and Experimental Results
The interference spectra at different positions were obtained in the measurement process as shown in Figure 7. The data in the upper left corner are the distance between the measurement signal and the reference signal, and the data in the upper right corner are the measurement result obtained by the interferometer. The larger the interval between the two signals, the higher the carrier frequency of the interference signal and the denser the interference fringes.
Because the spectrometer can only collect the interference signal by equal wavelength interval sampling in the measurement process, if the traditional time-frequency analysis method is used in the demodulation process of the interference signal carrier frequency, it will bring theoretical error. Because the spectrometer can only collect the interference signal by equal wavelength interval sampling in the measurement process, if the traditional time-frequency analysis method is used in the demodulation process of the interference signal carrier frequency, it will bring theoretical error.
Next, the error is analyzed simply according to the equation c = λf, where c is the speed of light, λ is the wavelength, and f is the optical frequency. Since the sampling interval of the spectrometer is equal to the sampling interval for the wavelength, and assuming the sampling interval of its wavelength is ∆λ, the corresponding frequency sampling interval can be expressed by the following equation: According to the Equation (26), when the spectrometer meets the equal wavelength interval sampling, the frequency interval sampling is not even. Although the traditional Fourier transform analysis method works well in processing continuous and uniform sampling time series, so that we could obtain a perfect spectrum, when dealing with the actual observation data that are non-uniformly sampled and contain a lot of noise, due to the non-uniformity of the data interval and finite time span, the power spectrum of Fourier transform would contain noise too, which would lead to the reduction of the characteristic spectrum power and the generation of a false spectral peak value of the signal. In addition, due to the enhancement of the noise signal, the amplitude and phase of the real periodic signal will be greatly affected, resulting in certain deviation, and in serious cases it will lead to failure of the parameter extraction. At the same time, to overcome the limitations in resolution, we adapted the Lomb-Scargle algorithm, which is based on the discrete Fourier transform principle: the time sequence is decomposed into a linear combination of some sinusoidal functions Y = A cos(ωt) + B sin (ωt), and then, the characteristics of the signal in the time domain can be transformed into the frequency domain.
In the Lomb-Scargle algorithm, the model curve of the data is fitted by using the sine function and least square method and the root mean square error (RMSE) is used to judge the coincidence degree of the implied periodicity trend of the data and the conjecture Next, the error is analyzed simply according to the equation c = λf, where c is the speed of light, λ is the wavelength, and f is the optical frequency. Since the sampling interval of the spectrometer is equal to the sampling interval for the wavelength, and assuming the sampling interval of its wavelength is ∆λ, the corresponding frequency sampling interval can be expressed by the following equation: According to the Equation (26), when the spectrometer meets the equal wavelength interval sampling, the frequency interval sampling is not even. Although the traditional Fourier transform analysis method works well in processing continuous and uniform sampling time series, so that we could obtain a perfect spectrum, when dealing with the actual observation data that are non-uniformly sampled and contain a lot of noise, due to the non-uniformity of the data interval and finite time span, the power spectrum of Fourier transform would contain noise too, which would lead to the reduction of the characteristic spectrum power and the generation of a false spectral peak value of the signal. In addition, due to the enhancement of the noise signal, the amplitude and phase of the real periodic signal will be greatly affected, resulting in certain deviation, and in serious cases it will lead to failure of the parameter extraction. At the same time, to overcome the limitations in resolution, we adapted the Lomb-Scargle algorithm, which is based on the discrete Fourier transform principle: the time sequence is decomposed into a linear combination of some sinusoidal functions Y = A cos(ωt) + B sin(ωt), and then, the characteristics of the signal in the time domain can be transformed into the frequency domain.
In the Lomb-Scargle algorithm, the model curve of the data is fitted by using the sine function and least square method and the root mean square error (RMSE) is used to judge the coincidence degree of the implied periodicity trend of the data and the conjecture model [18][19][20][21][22]. Consequently, when we use the Lomb-Scargle algorithm, Fourier transform can be applied to non-uniform sampled signals equivalently, which can not only contribute to extract weak periodic signals from the time series, but also to reduce the generation of false signals of the non-uniform time series to some extent. Therefore, the Lomb-Scargle algorithm was selected as the solution method in this experiment. The Lomb-Scargle algorithm is specially used for processing non-uniform sampling signals, and the results acquired also have better unimodality. Moreover, it has been proven that the Lomb-Scargle algorithm can demodulate the non-uniform sampling signal. The algorithm can be described by the following equations: where x j is the jth data, t j is the jth time, τ is specified for each ω to ensure time-shift invariance, and P(ω) is the power spectral density of the measured data we obtain with the unequal sampling interval for the wavelength. Figure 8 shows the comparison between the frequency spectrum of the interference signal obtained by Fourier transform and the spectrum obtained by the Lomb-Scargle algorithm. Compared with Fourier transform, the result obtained by the Lomb-Scargle algorithm has a better unimodal property. At the same time, compared with the measurement results of the interferometer, this data processing method has high precision.
Sensors 2022, 22, x FOR PEER REVIEW 10 of 12 model [18][19][20][21][22]. Consequently, when we use the Lomb-Scargle algorithm, Fourier transform can be applied to non-uniform sampled signals equivalently, which can not only contribute to extract weak periodic signals from the time series, but also to reduce the generation of false signals of the non-uniform time series to some extent. Therefore, the Lomb-Scargle algorithm was selected as the solution method in this experiment. The Lomb-Scargle algorithm is specially used for processing non-uniform sampling signals, and the results acquired also have better unimodality. Moreover, it has been proven that the Lomb-Scargle algorithm can demodulate the non-uniform sampling signal. The algorithm can be described by the following equations: where xj is the jth data, tj is the jth time, τ is specified for each ω to ensure time-shift invariance, and P(ω) is the power spectral density of the measured data we obtain with the unequal sampling interval for the wavelength. Figure 8 shows the comparison between the frequency spectrum of the interference signal obtained by Fourier transform and the spectrum obtained by the Lomb-Scargle algorithm. Compared with Fourier transform, the result obtained by the Lomb-Scargle algorithm has a better unimodal property. At the same time, compared with the measurement results of the interferometer, this data processing method has high precision.  As shown in Figure 9, different measuring points were selected in the range of 70 m. The experimental results show that the absolute distance measurement can be realized by dispersive interferometry in the range of 70 m, and the error of measurement was less than 12 µm. As shown in Figure 9, different measuring points were selected in the range of 70 m. The experimental results show that the absolute distance measurement can be realized by dispersive interferometry in the range of 70 m, and the error of measurement was less than 12 μm.

Conclusions
In this paper, we firstly analyzed different mathematical models of optical frequency comb with an odd or even number of teeth, and derived a unified mathematical model of the optical frequency comb signal, which is not affected by the number of comb teeth. Moreover, we derived the principle of absolute distance measurement based on DPI and proved the feasibility theoretically. At the same time, we compared the frequency spectrum of the interference signal obtained by Fourier transform and the spectrum obtained by the Lomb-Scargle algorithm, and then the latter was used to reduce the error caused by the non-uniform sampling of the spectrometer, which was verified by experiments.
A long absolute distance measurement experiment with mode-locked OFC was carried out; to overcome the limitations in resolution, the Lomb-Scargle algorithm was adopted to process the data. The distance resolution and the non-ambiguity range of the system are about 10 μm and 0.6 m. The experimental results show that the measurement error was less than 12 μm at the selected measurement points. The results show that the Lomb-Scargle algorithm can effectively reduce the error caused by the resolution, and the error of absolute distance measurement was less than 12 μm in the distance of 70 m based on the mode-locked OFC. Furthermore, owing to the advances in micro-resonator fabrication, there has been a great development in the research of chip-scale soliton microcombs (SMCs) with high repetition frequency, and in the future, we could combine the SMC with the Lomb-Scargle algorithm to eliminate dead zones based on DPI, which will be extended to arbitrary range measurement, as well as applied in ultra-precision semiconductor manufacturing, large aircraft processing, and other precision measurement fields.

Conclusions
In this paper, we firstly analyzed different mathematical models of optical frequency comb with an odd or even number of teeth, and derived a unified mathematical model of the optical frequency comb signal, which is not affected by the number of comb teeth. Moreover, we derived the principle of absolute distance measurement based on DPI and proved the feasibility theoretically. At the same time, we compared the frequency spectrum of the interference signal obtained by Fourier transform and the spectrum obtained by the Lomb-Scargle algorithm, and then the latter was used to reduce the error caused by the non-uniform sampling of the spectrometer, which was verified by experiments.
A long absolute distance measurement experiment with mode-locked OFC was carried out; to overcome the limitations in resolution, the Lomb-Scargle algorithm was adopted to process the data. The distance resolution and the non-ambiguity range of the system are about 10 µm and 0.6 m. The experimental results show that the measurement error was less than 12 µm at the selected measurement points. The results show that the Lomb-Scargle algorithm can effectively reduce the error caused by the resolution, and the error of absolute distance measurement was less than 12 µm in the distance of 70 m based on the modelocked OFC. Furthermore, owing to the advances in micro-resonator fabrication, there has been a great development in the research of chip-scale soliton microcombs (SMCs) with high repetition frequency, and in the future, we could combine the SMC with the Lomb-Scargle algorithm to eliminate dead zones based on DPI, which will be extended to arbitrary range measurement, as well as applied in ultra-precision semiconductor manufacturing, large aircraft processing, and other precision measurement fields.